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Abstract 

We comment on the results of a complete leading-twist next-to-leading or- 
der QCD analysis of the spacelike pion electromagnetic form factor at large- 
momentum transfer Q. For the asymptotic distribution amplitude, we have 
examined the sensitivity of the predictions to the choice of the renormalization 
scale. The results show that, regarding the size of the radiative corrections, 
reliable perturbative predictions for the pion electromagnetic form factor can 
already be made at a momentum transfer Q of the order of 5 to 10 GeV. 
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We comment on the results of a complete leading-twist next-to-leading order QCD 
analysis of the spacelike pion electromagnetic form factor at large-momentum 
transfer Q. For the asymptotic distribution amplitude, we have examined the 
sensitivity of the predictions to the choice of the renormalization scale. The re- 
sults show that, regarding the size of the radiative corrections, reliable perturbative 
predictions for the pion electromagnetic form factor can already be made at a mo- 
mentum transfer Q of the order of 5 to 10 GeV. 



1 Introduction 

The study of exclusive processes at large-momentum transfer represents a 
challenging area for application of perturbative QCD (PQCD). Although the 
PQCD approachtl undoubtedly represents an adequate and efficient tool for 
analyzing exclusive processes at very large momentum transfer, its applica- 
bility to these processes at experimentally accessible momentum transfer has 
been long debated and attracted much attention. In a moderate energy re- 
gion (a few GeV) soft contributions could still be substantial! Further on, 
the self-consistency of the PQCD approach was questioned regarding the non- 
factorizing end-point contributionsB It has been shown, that the incorpora- 
tion of the SudfJjCQv suppression in the so-called modified hard-scattering ap- 
proach (mHS A) &EI effectively eliminates soft contributions from the end-point 
regions and that the PQCD approach to the pion form factor begins to be self- 
consistent for a momentum transfer of about > 4 GeV^. However, in the 
PQCD approach to exclusive processes one still has to check its self-consistency 
by studying radiative corrections. It is well known that, unlike in QED, the 
leading-order (LO) predictions in PQCD do not have much predictive power, 
and that higher-order corrections are important. They have a stabilizing effect 
reducing the dependence of the predictions on the schemes and scales. 

In our recent jpapers □ we have clarified some discrepancies between pre- 
vious calculationsP and by including the complete closed form for the NLO 
evolution of the pion distribution amplitude (DA) derived recently^ we have 
obtained the complete NLO PQCD prediction for the pion electromagnetic 
form factor (within the so-called standard hard-scattering approach (sHSA) 
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and using different candidate DAs). The size of the NLO correction as weh 
as the size of the expansion parameter, i.e. QCD running couphng constant, 
serve as sensible measures of the self-consistency of PQCD prediction. But, 
as the truncation of the perturbative expansion for the pion form factor at 
finite order causes the residual dependence of the prediction on the choice of 
the renormalization and factorization scales (as well as on the renormaliza- 
tion scheme), the choices for these scales represent the major ambiguity in the 
interpretation of the results. 

In this paper we would like to outline our calculation^ and to address the 
scale ambiguity problem. 

2 Pion electromagnetic form factor in the sHSA 

In leading twist, the pion electromagnetic form factor can be expressed by a 
convolution formula 

FAQ^)^ f dx f dy <l>*{y,fij,)TH{x,y,Q^,fil,fil) <l>ix,fil). (1) 
Jo Jo 

Here = —q^ is the momentum transfer in the process and is supposed 
to be large, /z/j is the renormalization scale, and iip is the factorization scale 
at which soft and hard physics factorize; x and y {x = 1 — x and y = 1 — y) 
denote incoming and outgoing quark (antiquark) momentum fractions. 

The (process dependent) hard-scattering amplitude Th{x, y, Q^, /i^^) is 
calculated in perturbation theory and represented as a series in the QCD run- 
ning coupling constant asifJ-ji)- have used the dimensional Jjegularization 
method and the MS renormalization scheme in our calculation^ 

The (process independent) pion DA (p is intrinsically nonperturbative 
quantity, whose evolution can be calculated perturbatively and .represented 
as a series in as(/i|n). There are compelling theoretical results a^tl which dis- 
favor the end-point concentrated distributions, and one expects that the pion 
DA does not differ much from the asymptotic form. In this work we comment 
only on the results obtained with the asymptotic distribution (f>asix, ^J-'p)- 

Generally, one can express the NLO form factor as 

FAQ': mf) = Fi"^ (g2, f,l, ^,%) + (g2, ^% ^|,) . (2) 

The first term in (||) is the LO contribution, while the second term is the NLO 
contribution coming from the NLO correction to the hard-scattering amplitude 
as well as arising from the inclusion of the NLO evolution of the DA. For the 
results obtained using (pas{x, fijp) distribution, the effect of the NLO evolution 
of the DA is neghgible (« 1%). 
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Truncation of the pcrturbative series of Ft^{Q^) at any finite order causes a 
residual dependence on the scheme as well as on the scales (which is explicitly 
denoted in Eq. (||)). As we approximate only by two terms of the 

perturbative series, we hope that we can minimize higher-order corrections by 
a suitable choice of and /xf, so that the LO term -Fi"^ (Q^, Mfl; A^f) gives a 
good approximation to the complete sum 

We take that a PQCD prediction for pion form factor can be considered 
reliable provided the corrections to the LO prediction are reasonably small 
(< 30%) and the expansion parameter (effective QCD coupHng constant) is 
acceptably small (as(M|;) < 0-3 or 0.5). The consistency with the experimental 
data is not of much use here since reliable experimental dataE3 for the pion form 
factor exist for < 4 GeV^ i.e., outside the region in which the perturbative 
treatment based on Eq. (0)iii justified. It should also be mentioned that there 
are controversial argumentsEil regarding the reliability of existing experimental 
datatj The new data in this energy region are expected from the CEBAF 
experiment E-93-021. 

3 Examining the scale dependence of the NLO prediction 

The simplest and widely used choice for the and ^.p scales is 

^il^^Jil^Q\ (3) 

the justification for the use of which is mainly pragmatic. The prediction for 
the pion form factor depends very weakly on the choice of the factorization 
scale jip- Actually, taking to be an effective constant, i.e., = (Mf)j 
the only dependence of the results obtained using (pasix^fJ.'^) distribution 
comes from the NLO evolution of the DA which is negligible. 

Physically, a more appropriate choice for would be that corresponding 
to the characteristic virtualities of the particles in the parton subprocess (which 
are considerably lower than the overall momentum transfer i.e., virtuality 
of the probing photon) 

nl^a{x,y)Q\ (4) 
For example, some of the physically motivated choicesi are 

a{x,y) e {xy ,^xyy .e'^/'^xy] . (5) 

These correspond, respectively, to the (LO) gluon virtuality, geometrical mean 
of the gluon and quark virtualities (an attempt to take into account that the 
QCD coupling is renormalized not only by the vector particle propagator, but 
also by the quark-gluon vertex and the quark-propagator), and to the choice of 
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X=<l-x>=l/2 



-— Mr = X q 
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Figure 1: NLO prediction for Ftt(Q^) obtained using the 4>aa{x, y?p) distribution ampli- 
tude and the choices of given by Eqs. (^), and (^), while = and (x) = 0.5. 
The shaded area denotes the theoretical uncertainty introduced by the renormalization scale 
ambiguity. 



the renormalization scale according to the Brodsky-Lepage-Mackenzie (BLM) 
procedure E3 the essence of which is that all vacuum-polarization effects from 
the QCD (3 function should be resummed into the running coupling constant 
(the NLO coefficient of becomes n/ (i.e. (3q) independent). 

A glance at Eq. (|i|), where the coupling constant as(p\) appears under 
the integral sign, reveals that the choice (||) leads immediately to the problem 
if the usual one-loop formula for the effective QCD running coupling constant 
is employed. Namely, regardless of how large is, the integrations over x and 
y allow oisit^'R) to be evaluated at low momenta i.e., in the region where usual 
one-loop formula is not a good recp^sentation of the effective QCD coupling. 
There are number of proposals PrM for the form of the coupling constant 
Q!5(/i|j) for small but its implementation in this calculation demands the 
more refined treatment^ Alternatively, one can choose to be an effective 
constant 

MH = (Mfl) = («(a^,y)) Q' = aQ'- (6) 
Hence, the expressions (||) get replaced by their respective averages 



(7) 



Now, the key quantity in the above expressions is (x ) , the average value of the 
momentum fraction. Owing to the fact that (f>as{x, is centered around the 
value X = 0.5, the simplest choice is {x) = 0.5. 
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Figure 2: The ratio Fi^\Q^)/F^\Q^) obtained using the same DA and the same choices 
for fj,^ and fi^ scales as in Fig. hi The shaded area denotes the region of predictions which 

corresponds to l^i^^ (Q^)/^^"' (Q^)! < 30%, and as(/.i|j) <0.5 (as (mIj) < 0.3). 

Numerical results of our complete NLO QCD calculation for Q^Ft^{Q^), ob- 
tained using the (j)as{x, fJ,%) distribution, with ul^ = and different choices for 
the renormalization scale /ifj given by (j^) (||), (0), and {x) = 1/2, are displayed 
in Fig. |l| (in our calculation we take Ajjg — 0.2 and f-^ = 0.131 GeV). The 

ratio of the NLO to the LO contribution to F^(Q^), i.e., Fi^\Q^)/ Fi°\Q^), 
as a useful measure of the importance of the NLO corrections, is plotted as a 
function of in Fig. ^. 

One notices that the total NLO perturbative prediction for Q'^Ft^{Q'^) is 
somewhat below the trend indicated by the presently available experimental 
data, but what alarms us and seems to question the self-consistency of the 
PQCD approach is the fact that the ratio is^^-'/i^r^-' corresponding to the 
often encountered choice /i^ = (solid line) is rather high: F'ip / F^^ < 30% 
is not reached until « 500 GeV^. The answer to this problem lies in the 
previously stated inappropriateness of the choice = Q^. Namely, owing 
to the partitioning of the overall momentum transfer among the particles 
in the parton subprocess, the essential virtualities of the particles are smaller 
than Q^, so that the "physical" renormalization scale, better suited for the 
process of interest, is inevitably lower than Q^. It follows from Fig. ^ that 
by choosing the renormalization scale determined by the dynamics of the pion 
rescattering process, the size of the NLO corrections is significantly reduced and 
reliable predictions are obtained at considerably lower values of Q^, namely, 
for Q2 < 100 GeV2. 



' ^^^^ X=<l-x>=l/2 


^^ = Q 

^^ = X^Q' 

— - m^ = x'"q^ 
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The total NLO prediction for Q^F.^{Q^) depends weakly on the choice of 
fij^. This is a reflection of the stabilizing effect that the inclusion of the NLO 
corrections has on the LO predictions. Let us explore this point more closely. 
By taking and /J^ to be effective constants, the LO and NLO contributions 
to the pion form factor obtained using (f>asix, fi^) distribution amount to 



Q^F^°\Q',^,l,^^%)=8 7Tfias{^^l), (8) 



(9) 



where the NLO contribution coming from the NLO evolution of the DA is ne- 
glected. As it is seen from (||), all of the /iij dependence of the LO result Q'^f!^^ 
is contained in the strong coupling constant as{p\)- Thus, as decreases the 
LO result increases, and it increases without bound. In contrast to the LO, the 
NLO contribution Q'^Fi^\ as evident from the explicit expression given in (H), 
decreases (becomes more negative) with decreasing ^lj^. Upon adding up the 
LO and NLO contributions, we find that the full NLO result, as a function of 
stabilizes and reaches a maximum value for — ^i'ixtreme ~ Q^/IS. If we 
take that the renormalization scale continuously changes in the interval defined 
by a G [1/18, 1] the curves representing the NLO predictions for Q^FTr{Q^) fill 
out the shaded region in Fig. |l], and this shaded region essentially determines 
the scale ambiguity related theoretical uncertainty of the NLO calculation. 



4 Towards resolving the renormalization scale ambiguity problem 

The optimization of the scale and scheme choice according to some sensible 
criteria remains an important task for the application of PQCD. Several scale- 
setting procedures were pr opo sed in the literature: the principle of fastest 
apparent convergeare (FAC)lla (Fi^' [Q^, /^fac) ~ 0), the principle of minimal 
sensitivity (PMS)ll3 (/ipMS = l^extreme)^ the BLM method. The applica^ 
tion of those methods can give strikingly different results in some calculationstH 
As it is known, the relations between physical observables must be inde- 
pendent of renormalization scale ancLscheme conventions to any fixed order of 
perturbation theory. It was argued Efl that applying the BLM scale-fixing to 
perturbative predictions of two observables in, for example, MS scheme and 
then algebraically eliminating ajjg one can relate any perturbatively calcula- 
ble observables without scale and scheme ambiguity, where the choice of BLM 
scale ensures that the resulting "commensurate scale relation" (CSR) is inde- 
pendent of the choice of the intermediate-renormalization scheme. Following 
this approach, in paper by Brodsky et aliB the exclusive hadronic amplitudes 
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Table 1: NLO PQCD results for the pion form factor, F-^{Q^) in the ay scheme. 



[GeV'] 




Fy^(Q2)/Fr(g2) [%] 


Q^F^Q^) [GeV^] 


20 
100 


0.434 
0.289 


-26.5 
-17.7 


0.138 
0.103 



were calculated in ay scheme, in which the effective coupUng Q!y(/i^) is defined 
from the heavy-quark potential V{n^). 

It follows from that the Po dependent term vanishes for fij^ — {lJ.%i^j^.f) = 
e^^^^'^Q^. Using the scale-fixed relationEj between ay(^y) and ajjg{fj^^) = 
asipji), and Eqs. (SP) one obtains the NLO prediction for the pion form 
factor in ay schemeEj 

g2^^.(Q^M?^) = 8^/>v(M^v')(l-l•92^^) , (10) 

where = e^^^fi\ = e^^Q"^ k, Q-^/20. Considering the energy region we 
are interested in, for the purpose of this calculation we approximate ay with 
usual one-loop expression. Thus obtained numerical predictions are given in 
Table |l|. Considering the size of the i^i^V^i°^ ratio (« -20%) and the size 
of the effective coupling, we find that NLO PQCD predictions for the pion 
form factor obtained in ay scheme can be considered reliable for below 100 
GeV^ i.e., already for momentum transfer Q of the order of 5 - 9 GeV. 



5 Conclusions 

We have shown that, regarding the size of the radiative corrections, the sHSA 
can be consistently applied to the calculation of the pion electromagnetic form 
factor already for momentum transfer Q of the order of 5 - 10 GeV. In order 
to improve this PQCD based prediction it is necessary to obtain and apply 
the proper form of the QCD coupling in the low-momentum regime as well as 
to investigate the corrections (present in the few GeV region) introduced by 
the mHSAtj Reliable experimental data and inclusion of the soft contributions 
should then enable complete confrontation between theory and experiment. 
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